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STUDY BY THE PRANDTL- GLAUERT METHOD OF COMPRESSIBILITY 
EFFECTS AND CRITICAL MACH NUMBER FOR ELLIPSOIDS OF 
VARIOUS ASPECT RATIOS AND THICKNESS RATIOS 
By Robert V. Hess and Clifford S. Gardner 


SUMMARY 


By the use of a form of the Prandtl-Glauert method 
that is valid for three-dimensional flow problems, the 
value of the maximum incremental velocity for compressible 
flow about thin ellipsoids at zero angle of attack is 
calculated as a function of the Mach number for various 
aspect ratios and thickness ratios. The critical Mach 
numbers of the various ellipsoids are also determined. 

The results indicate an increase in critical Mach number 
with decrease in aspect ratio which is large enough to 
explain experimental results on low- aspect-ratio wings 
at zero lift. 


INTRODUCTION 


Recent tests (references 1 and 2) have shown that 
an appreciable increase in the critical Mach number, 
together with other improvements of the aerodynamic 
characteristics at supercritical- Mach numbers, results 
from the use of wings of very low aspect ratio. These 
improved characteristics have been somewhat qualitatively 
ascribed to "three-dimensional relief," although no 
quantitative theoretical discussion has yet been provided. 

In the present paper an effort is made to provide 
such a study by considering the flow, at zero angle of 
attack, about a series of thin ellipsoids of various 
aspect ratios and thickness ratios. Ellipsoids were 
chosen because they are amenable to calculation. 
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Although they differ appreciably from the wings of 
reference 1, which had an NACA 0012 airfoil section and 
rectangular plan form, ellipsoids should nevertheless 
show similar aspect-ratio effects. The calculations 
were made for ellipsoids of thickness ratios 0.10, 0.15, 
and- 0.20, and for the entire range of aspect ratios 
from the elliptic cylinder to the ellipsoid of revolution. 

The compressibility effects were computed by the use 
of a form of the Prandtl-Glauert method that is valid for 
three-dimensional flow problems. The method has been 
given by Go the rt (reference 5)' without, however, very clear 
mathematical proof. Since the methods that have been 
commonly used (see, for example, reference Ip, 5, and 6) 
are applicable only to two-dimensional problems, a detailed 
proof of the method correct for three-dimensional flow is 
included in the - appendix,. A brief discussion of 1 he 
accuracy of the Prandtl-Glauert method, as applied to 
ellipsoids, is also given. 

This study was made during the period from October 
191+5 to April 19I1.6. 


SYMBOLS 

U free-stream velocity 


C velocity of sound in free stream 

M free-stream Mach number (U/C) 

Y ratio of specific heats 


P = Vl - M 2 

x, y, z rectangular coordinates 


B thin body 

cp velocity potential 


u, v, w x-, y- , and z -components of incremental 
velocity for compressible flow about B 


B» 


body obtained by stretching B in direction 
of x-axis by the factor l/S 


u 1 , v', w 1 x-, y-, and z-components of incremental 

velocity, for incompressible flow about B' 
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a 

b 

c 

a* - 


a 


u - Umax 
U 


u(M) 

value of 
to M 

u(0) 

value of 

€ 

thicknes 


maximum semichord of ellipsoid 

semispan of ellipsoid 

maximum semi thickness of ellipsoid V- 


a. = - ^ - ^ -1 — = k 

rrab n a / 




value of u for incompressible flow (M = 0) 
... ( Thickness 


■) 


jj\ ^ y Hi / 


§U,o) 


value of ratio of incremental velocity to 
free.- stream velocity for compressible 
flow having Mach number M about a body 
having thickness ratio € 

value of ratio of incremental velocity to 
free-stream velocity for incompressible 
flow about a body having thickness ratio € 


Subscript: 

max maximum value 


METHODS 0? CALCULATION 


The Pr andtl -Glauer t method for three-dimensional 
flow . - The Prandtl-Glauer t method is used in the present 
paper in the following form: The incremental velocities 

at a point P on the surface of a thin body B in three- 
dimensional compressible flow ma? be obtained, in three 
steps: 

(l) The x-coordinates of all points of B are 

increased by the factor l/p, where $ = yfl - ¥?- and 
the x-axis is in the stream direction. This trans- 
formation takes B into a ‘'stretched" body B 1 . 
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(2) The velocity increments u', v’, and w’ 
parallel to the x-, y-, and z-ax.es, respectively, at 
the point. P' on the stretched body B' corresponding 
to the point P on the original body B are calculated 
as though B ’ were in an incompressible flow having the 
same free-stream velocity as the original compressible 


flow. 


( 3 ) The values u, v, and w of the incremental 
velocities at the point P on B in compressible flow 
are then given by the equations 

u = u * 

P 2 

v = i v ' 

• P 

• w = — w ’ 

P 

A derivation of this form of the Prandtl-Glauert 
method is given in the appendix. The method in essentially 
this form has been given by Gother t (reference 3 ) without, 
however, a very clear proof. ( Gother t prefers to shrink 
the lateral coordinates of the body by the factor (3 
rather than to expand the coordinate in the stream 
direction by the factor l/p; obviously the two pro- 
cedures lead to the same result. ) Prandtl (reference I4.) 
and von Harman (reference 5) state the method in a form 
that is valid for two-dimensional flows but in general is 
incorrect for three-dimensional flows. Goldstein and. 

Young (reference 6) also give a discussion leading to 
results that are correct only for two dimensions. A 
discussion of the reasons for the failure of these ' 
commonly used methods for three-dimensional flow problems 
is included, in the appendix. 


Calculation of incremental velocity for compressible 
flow about ellipsoids.- In order to determine, by the 

the incremental velocity on the 


a , b , and c 


Prandtl-Glauert method, 
surface of an ellipsoid having semiaxes 
where a is the length of the semi axis in the stream 
direction, the incremental velocity is calculated for 
stretched ellipsoid having semiaxes a', b, and. 


a 


c , where 


a’ = •£, in an incompressible flow having the same 
stream velocity, and the result is multiplied by 1 /( 3 ^ 
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For incompressible ; f' : l'ow about the stretched ellipsoid,' the 
velocity potential on the surface of the ellipsoid is 
given by 


9 - TJx 

2j - a r, 


where 


a 0 = a’ be 


d/\ 


v 0 


( a 1 2 + X ) _v/( a ’ ^ + a . ) ( + X. ) ( c^ + X. ) 


(see, for example, reference 7 ) • The incremental velocity 
at x ="0 (half -Chord line on the -stretched ellipsoid in 
incompressible flow) is then sriven bp 

-- • -j 


• _ ■ : ■. • 2. - a Q . 7/ .... 

This value is .the maximum value of u* (reference 7) 
and evidently is the same at all points on the half-chord 
line. The incremental velocity at the half -chord line for 
the compressible flow about the original Unstretched ellipsoid 

is given by 


u 




U 


(1) 


Various formulas are necessary for the evaluation of the 
integral a Q when a’ > b > c, b > a* >c, or a' > b = c 
(ellipsoid, of revolution). 

For a' > b > c, the value of a Q is given by the 
formula 

• . . ■ ■ a 0 - - -- -- - (F - e) (2) 

(a'2 - b^) \/a ,<;1 - c^ ' 

where F and E are incomplete elliptic integrals of the 
first and second kind, respectively, defined as follows: 


F 



/ 

Mo 
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where 


COINCIDENT! Aft NACA RK No. LjBOpa 


E 



sin^\}/ d\{/ 


k 


, /a ’ 2 - b 2 
\J a’ ^ - c 2 


and 


sin cp 


_ Va 1 2 c 2 


' For b > a* •> c, the value of a Q is • given by . the 
formula 


2 a , bc\/b^ - c 


a - — — — ; — i n 

. (h? - a' 2 ) (a '2 - = 2 ) j. 


/• 


v b 2 - c 2 


a 


-i 2 -* (5) 


where F And E are defined as before,, with 


k = 


■V 


fb 2 - a ' 2 


b 2 - c 2 


and • 

/ - “p' 

sin cp = V'k~ ., 

. b 

Equation (2) is derived .from the first equation given 
in equations ( 5 . 1 , 3 ) of reference 7 by substituting a' for 
a and by using the expression for. k in terms of a', b, 
and c. Equation (3) is derived from the second equation 
given in equations ( 5 . 13 ) of reference 7 hy^ interchanging a 
and b, substituting a' for a, • and using the expression 
for k in tepms of a', b, and c. 

For a' > b = c, (ellipsoid of revolution), a Q is 
given by the equation 

' p°°' 

_ , , p ■ • ci ^ 

°° = y a .2 + X >/ 2 ( b 2 +x ) 
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which resolves into 

a o = ~ (log e -■ — -• 2e) 

e 5 V 1 - e / 

where 


e = - b2 

a* 


If this value' for a D is. substituted in equation (1), 
the incremental velocity at the half-chord line for the 
ellipsoid of revolution is found to be 


u 


1 

P" 


log e - ■ t " e - 2e 

l - e 


2e 


1 - e 


2 


l°g e 


1 -f G 
1 - ^ 


The limiting case cf infinite aspect ratio (elliptic 
cylinder) was treated by the use' of formulas for' the 
ellipse in two-dimensional flow (•reference 8). 

Calculation of the critical Mach number . - For flow 
about a two-dimensional body" the fres-stream Mach number 
for which sonic speed is first reached at some point on 
the surface is called the critical Mach number / because 
of the development of shocks and the accompanying deteri- 
oration of the aerodynamic characteristics shortly after 
this Mach number is exceeded. For the general three- 
dimensional body, however, exceeding the Mach number at 
which sonic speeds first appear does not necessarily 
imply the possibility of shock formation; for example, 
in the case of the infinite yawed cylinder (reference 9) 
shocks may be impossible even when supersonic speeds 
exist cn the si?rface. In this case, since the flow 
must be the same at corresponding points of the sections 
along the cylinder, any shock front must be parallel to 
the axis of the cylinder. Such an oblique shock requires 
that the velocity component normal to the shock front - 
that- is, normal to the cylinder axis - be supersonic. 

The magnitude of the velocity component parallel to the 
axis is Immaterial. 

By what seems a reasonable extension. of this concept, 
the following criterion for the possibility of shock 
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formation on the general three-dimensional body is 
tentatively proposed: Consider the line on the body 

that forms the locus of those points where the total 
velocity on each streamline is a maximum (that is, 
where the pressure is a minimum); a shock will form 
when the velocity normal to this line exceeds sonic 
velocity at some point along the line. 

It should be emphasized that, for the three- ( 
dimensional body, the shock may begin to form over 
only a very small region of the surface, so that, in 
general, existence of the condition just defined does - 
not necessarily imply an imminent deterioration of the 
aerodynamic characteristics of .the body. 


For the special case of the unyawed ellipsoid 
considered in the present paper, however, no appreci- 
able' analysis of shock formation or shock extent along 
the lines just indicated seems to be required. As. is 
shown in the section "Calculation of Incremental Velocity 
for Compressible Flow about Ellipsoids" , the maximum 
velocity for an unyawed ellipsoid is in the stream 
direction and occurs simultaneously at all points along 
the half-chord line. Sonic velocity is thus reached 
simultaneously along a line that extends across the 
entire span of the body and is normal to the stream 
direction. These conditions also exist in the case of 
the unyawed infinite cylinder, that is, the two-dimensiona 
body. 


The critical Mach number of the ellipsoid was 
accordingly determined by solving graphically the 
equation 


u(M) 


1 



+ 





Y + 1 
2 


-1 


where u(M) is the ratio of the incremental velocity 
at the half-chord line to the stream velocity at the" 
Mach number M. 


Accuracy of the Prandtl-Glauert method .- The 
Prandtl-Glauert method is based on the assumption of 
small perturbations. Consequently, near the" nose of 
the ellipsoids discussed in the present paper, where 
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the assumption of small perturbations is violated, the 
results given by the Prandtl-Glauer t method cannot be 
expected to be. reliable. More reliable values, however, 
should be obtained for the maximum incremental velocity, 
which occurs at the half-chord line. The accuracy of 
the Prandtl-Glauert approximation for the maximum incre- 
mental velocity may be estimated by cdmparison with more 
exact solutions of the compressible flow problem. An 
iteration method in which the f randtl-Glauert method is 
used as the first approximation has been proposed by 
Busemann (reference 10 ). The first and second approxi- 
mations have been calculated by Hantzsche and Wendt for 
the elliptic cylinder (reference 11) and by Schmieden 
and. Kawalki for the ellipsoid of revolution (reference 12). 
Calculation of the maximum incremental velocity for the 
elliptic cylinder having thickness ratio 0,20 by a 
formula for the second approximation given in reference 10 
shows that the value given by the Prandtl-Glauert method 
at a Mach number of 0.8 is almost 20 percent lower than 
the value given by the second approximation. For the 
ellipsoid of revolution, however, the va^ue of the 
maximum, incremental velocity given by the Prandtl- 
Glauert method agreed with the value given by the second 
approximation to within 5 percent at a Mach number of 0.8 
for thickness ratios up to C.JO. Although the second 
approximation is not the exact solution, ~it indicates 
that the error involved in using the Prandtl-Glauer t 
method to estimate the maximum incremental velocity for 
ellipsoids having a given thickness ratio is greatest 
for the limiting case, of -the elliptic cylinder (A = =») • 
and very small for the ellipsoid of- revolution, which 
has a very low aspect ratio. . The .error may be expected 
to be intermediate in magni tude. f or intermediate values 
of the aspect ratio and tc decrease with' aspect ratio. 

The reduction of error of the Frandtl— Glauert method with 
a decrease in aspect ratio was to be expected, as the 
incremental velocities are smaller for ellipsoids having 
low aspect ratio. 


RESULTS AND DISCUSSION 

Results . - Figures. 1, 2, and 3 show the value of the 
velocity ratio u = - at the half -chord line plotted 
against the Mach number for ellipsoids at zero angle of 
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attack for various aspect ratios and section thickness 
ratios equal to 0 . 10 , 0 . 15 , and 0 . 20 . In the same 
figures the sonic velocity boundary having the equation 



-1 


is plotted. The abscissa of the intersection of this 
boundary line with the curve of u plotted against M 
for any aspect ratio is the critical Mach number. In 
order to show the effect of compressibility more directly, 

the ratio of maximum, incremental velocity for com- 


pressible flow to the maximum incremental velocity for 
incompressible flow for the same free-strea.m velocity is 
plotted against the Mach number in figures ip, 5, and 6 
for the same aspect ratios and. thickness ratios. Similar 
curves for the ellipsoid of revolution, which is a 
special case of the ellipsoid, having three unequal axes, 
are plotted for the same thickness ratios in figures 1 
to 6. Figure 7 presents curves of critical Mach number 
against aspect ratio for thickness ratios of 0.10, 0.15, 
and 0.20. ' ’ 


Three-dimensional relief .- It may be seen from 
figures i. to 3 that the three-dimensional relief, that 
is, the difference between the velocity on the ellipsoid 
and the velocity on the corresponding ellipsoid of 
ini inite aspect ratio (elliptic cylinder), increases 
with a decrease in the aspect ratio. This increase has 
two causes: 

(1) For a^flow with M equal to zero (incompressible 
flow), the relief effect increases with a decrease in the 
aspect ratio. 

( 2 ) For larger values of M (compressible flow), an 
additional relief effect occurs with a decrease in the 
aspect ratio because of the fact that the compressibility 
effect (increase of incremental velocity with an increase 
in the Mach number) decreases with a decrease in the 
aspect ratio. (See figs. Ip to 6.) It may be seen that 
this additional three-dimensional relief increases most 
rapidly. at high Mach numbers. 
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Prom figures 1 .to 6 it may be seen that the compressi- 
bility effect on the maximum incremental velocity is 
greatest for A equal to infinity (infinite elliptic 
cylinder) and is smallest for the ellipsoid of revolution. 
The compressibility, effect on the maximum incremental 
velocity for the elliptic cylinder is proportional to 

1 

which is in agreement with the usual form of 

\jl - 

the Prandtl-Glauert method in two dimensions. The com- 
pressibility effect on the maximum incremental velocity 
for the ellipsoid of revolution is small in comparison 
with that of the elliptic cylinder. In fact, as the 
thickness ratio of any type of body of revolution 
approaches zero, the compressibility correction factor 
approaches unity, for in this limit the incremental 
velocity in incompressible flow is proportional to the 
square of the thickness ratio, so that the effect of 
stretching the body (first step of Prandtl-Glauert method, 
see the appendix) is exactly compensated for by the multi- 
plication of the incremental velocities by l/p2 (third 
step of the Prandtl-Glauert method). For ellipsoids of 
practical thickness ratios, however, the incremental 
velocity varies more slowly than the square of the thickness 
ratio. The compressibility effect for the ellipsoid of 
revolution (figs. 4, '5, ana 6 ) is thus considerable at 
high Mach numbers. For example, for a thickness ratio of 
0*20 and at a Mach number of 0,8, the compressibility effect 
amounts to about JO percent of the incremental velocity in 
incompressible flow. 

The effect of the thickness ratio on the three- 
dimensional relief may be seen by a comparison of figures 1, 
2, and J. From figure 1 it may be seen that, for a 
thickness ratio of 0.10, at a Mach number of 0.75, the 
maximum incremental velocity for A = 2 is 76 percent 
of the maximum incremental velocity for A = ». From 
figure on the other hand, it may be seen that, for 
a thickness ratio of 0.20, at a Mach number of 0.75, the 
maximum incremental velocity for A = 2 is 75 percent 
of the maximum incremental velocity for A = os. Thus, 
an increase in the thickness ratio causes only a very 
small increase in the three-dimensional relief. 

Critical Mach numbe r.- Figures 1, 2, J, and 7 
indicate chat ah increase in the critical Mach number 
of an ellipsoid at zero . Ilf t may be obtained by decreasing 

- OOMZCTTTIAb 
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the aspect ratio. For example, for ellipsoids having a 
tnickness ratio of 0.10, a decrease in the aspect ratio 
from infinity to 2 causes the critical Mach number to 
increase from 0 . 82 ? to O.857. For. a thickness ratio of 
C .20 , a decrease in the aspect ratio from infinity to 
2 causes the critical Mach number to increase from O.7I4.I 
to O.783. Although ellipsoids having greater thickness" 
ratio have lower critical Mach numbers, a decrease in 
the aspect ratio is slightly more effective in increasing 
critical Mach numbers for ellipsoids c-f greater 
thickness ratio. Figure 7 indicates that only a large 
reduction. in aspect ratio will cause a significant rise 
of the critical Mach number. 


Comparison with test re sults on low aspect ratio 
wings . - Figure 6 of reference 1 shows the minimum drag 
coefficient (C p for zero lift) plotted against the Mach 
number for. wings having an N AC A 0012 section and various 
aspect ratios. The critical Mach number for any asoect 
ratio may be estimated roughly as the Mach number for 
which the. drag coefficient first begins tc rise. The 
rough estimate of the critica.1 Mach numbers obtainable 
by this consideration is not. sufficiently accurate to 
warrant comparison of the numerical values with the 
numerical values of the critical Mach number obtained 
in the present pacer for thin ellipsoids. Comparison of 
the numerical results is, moreover, not warranted inasmuch 
as the wings of reference 1 did not have an elliptic 
section and furthermore had a rectangular plan form. 

A qualitative comparison may be made 7 however, between 
the results of the present paper and. those of reference 1. 
The increase in critical Mach number with decrease in 
aspect ratio indicated in figures 1,2, 3, and 7 of the 
present paper is considered sufficiently large to explain 
the corresponding effect indicated in figure 6, reference 1 


It is mentioned in reference 1 that the Mach number 
for a significant rise in the drag coefficient is approxi- 
mately 0.1 higher for an aspect ratio of 2 than for an 
infinite aspect ratio. This vslue is considerably higher 
than the increase. in critical Mach number due to a decrease 
m the aspect ratio. Since, for low-aspect-ratio win, os, 
the drag coefficient increases only gradually after the * 
critical Mach number is reached, the critical Mach number 
for a wing having low aspect ratio does not indicate so 
critical a change in the flow phenomena as the critical 
Mach number for a wing having high aspect ratio. It is 
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thought that the smaller rate of increase of the drag 
coefficient for wings having low aspect ratio is due to 
the fact that, at the critical Mach number, the rate of 
increase with Mach number of the incremental velocity is 
less than for high aspect ratios, as may be seen from 
figures 1, 2, and 3. 


CONCLUSIONS 


A study by the Prandtl-Glauert method of compressi- 
bility effects and critical Mach number for ellipsoids of 
various aspect ratios and thickness ratios indicated the 
following conclusions: 

1 • The flow about the unyawed ellipsoid is analogous 
to that- about the infinite unyawed cylinder in that 
sonic velocity is reached, simultaneously along a line 
that .extends across the entire span of the body and is 
normal to the stream direction. 

2. The critical Mach number for a thin ellipsoid 
may be predicted with good accuracy by means of the 
Prandtl-Glauert method, and the accuracy increases with 
decrease in aspect ratio. 

3. The compressibility effect on the flow about an 
ellipsoid decreases as the aspect ratio decreases. 

The three-dimensional relief for ellispoids is 
essentially independent of the thickness ratio, for 
thickness ratios from 0.10 to 0.20. 

5. For ellipsoids of thickness ratio 0.20, the 
critical Mach number increases by O.Oij. when the aspect 
ratio is changed from « to 2; for ellipsoids of thickness 
ratio 0.10 the increase is 0.03. 

6. The calculated increases in critical Mach number 
are sufficiently large to explain the experimentally 
observed increases in the Mach number at which the drag 
first begins to rise. 

7. The experimentally indicated reduced rate of 
drag rise for low-aspect-ratio wings at zero lift as 

, aamamm&sr 
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compared to that for wings having infinite aspect ratio 
may be explained qualitatively on the basis of the 
results obtained for the three -dimensional relief for 
ellipsoids. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee For Aeronautics 
Langley Field, Va. 
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APPENDIX 


THE PRANDTL-C-IAUERT METHOD FOR 


THREE - DIMENS ION AL FL OW 


.A derivation of the Prandtl-G-lauert method for three- 
dimensional flow .- A brief derivation of a. form of the 
Pr and 1 1 - G1 auer t method correct for three dimensions may • 
he given as follows: A first-order' approximation to the 

subsonic compressible flow about a thin body B, the 
surface of which has the equation 

S ( x , y , a ) = 0 


may be obtained by finding a solution of the linearized 
differential equation for the potential cp of the 
incremental velocities, 


«2 


CO 


XX 


( p. 


yy 



( A1 ) 


where the x-axis is in the stream direction and the 

incremental velocities cp„, © , and cp are small 

* y z 

compared with the stream velocity U. At all points on 
the surface of B, the potential cp must satisfy the 
boundary condition 


(U + 


<s> x ) s 


X 


Cp S + 
7 7 


Vz 


= 0 


(A2) 


which states that the flow is tangential to B. 
is assumed thin, ' S x is small compared with 3y 

consequently the second-order term cp S may be 

A A 

and the boundary condition becomes 


Since B 
and S z ; 

neglected. 


TJ£ 


9,. 


7 7 


+ cp. S = 0 

•Z Z 


In order to solve the boundary-value problem given 
by equations (Al) and (A2) in terms of incompressible 
flow the following transformation of variables is used: 



U3) 



l6 


C UiMi 1 ' i TK fT 


l 
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Under this transformation equations (Al) and (A 2 ) become, 
respectively. 


9 


1 

x ’x 1 


+ Cp 1 


yy 


9 



(All) 


US ,+cp’S + cp 1 S = 0 
x y y Z Z 


U5) 


Equations (Al+) jznd ( A5 ) are, respectively, the differ- 
ential equation and boundary condition for the potential 
9' of the incremental velocities of an incompressible 
flow with free-stream velocity U, in the x’ , y, 2 
space, about a thin body B' , the surface of which has 
the equation 


' S(px', y, z) =0 

The incremental velocities in the compressible flow 
are thus given by 


U = 9 Z = 

-i_ cp ' , = 

0 y x» 

h". 



p- 

p 


v = 9 = 

/ 

i Cp * = 1. 

P y p 

V * 


W = cp z = 

1 cp' = 1 

R 2 A 

W ' 



K P 

1 


where u, v, and w and 

u', v’, 

and w ' 

are the 


incremental velocities at corresponding points in the 
compressible flow about B and the incompressible flow’ 
about B', respectively. 


The foregoing analysis establishes the Prandtl- 
Glauert method for three-dimensional flow in the 
fol-lowing form: The incremental velocities at a point 

P on the surface cf a thin body B in compressible 
flow may be obtained in three steos: 

* y / 

( 1 ) The x-coordinates of all points of B are 
increased by the factor l/{ 3 , where 

P = Vi - M 2 

and where the x-axis is in the stream direction. This 
transformation takes B ..into a stretched -body B’. 
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(2) The incremental velocities u*, v ' , w’, in the 
direction of the x-, y-, and z-axes, respectively, at the 
point P 1 on 3' corresponding to the point P on 3 
are calculated as' though B* were in an incompressible 
flow having the same free-stream velocity as the original 
compressible flow. 


( 5 ) Tne values u, v, and. w of the incremental 
velocities at the point P on the original unstretched ■ 
body B in compressible flow are then "found bv the 
equations 



I 


Failure for three-dimensional flow pr oblems of the 
commonly stated for ms of the .Frsn.dti.-Gl anert. ' — 
According to the form of the Prandtl-Glauert method given 
by Prandtl (reference Ip) and von Karmen (reference 5 ), 
the incremental velocities for a compressible flow about 
a thin body E are the same as the incremental velocities 
ox corresponding points for incompressible flow having the 
same fro-e-stream velocity about a body obtained by expanding 


T? 


in the directions normal to the free-stream direction 


. w v — owj . O' cnu j. jl tr; u u J. LU1 

by the f ac uor l/p. That is, - for bodies of revolution, cr 
two-dimensional bodies, 


H ( e >.«) - u vb n \ 
•U ' ' u » °' 

According to. Gather t ’ s method, however, 


# (e, M) = ~ ~ (pe,‘ 0) 


( a 6 ) 


Thus, Prandtl 's and von Karma n 1 s method is valid only if 
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that is, if and only if the incremental velocity for 
incompressible flow about the bodies under consideration 
is oroportional to the thickness ratio. This relation 
is approximately valid for thin two-dimensional bodies, 
so that the method of Prandtl and von Karman may be 
expected to be 'valid for two-dimensional flows. The 
relation is not true in general for three-dimensional 
bodies; for example, for a very- thin body of revolution 
the incremental velocity is more nearly proportional to 
the square of the thickness ratio than to the first power 

/ 

Von Karman approaches the problem by making the 
transformation 

y' = Py 

z 1 = pz 

<p ' = cp 


Under this transformation the linearized equation of 
compressible flow goes into Laplace ' s equation; however, 
the transformed boundary condition is not satisfied c-n 
the surface of the transformed (contracted) body but on 
the surface of an expanded body. Thus, the boundary 
condition is not satisfied cn the boundary but at points 
near the boundary. This procedure is applicable to two- 
dimensional problems (as, for example, in the thin-wing 
theory, reference 1 5 ) , because the velocity increments 
induced by the equivalent line distribution of singu- 
larities vary only slowly in the neighborhood of the 
line of singularities. For a body of revolution, how- 
ever, the velocity increments induced by a line of 
singularities go to infinity at the line of singularities 
for such bodies, accordingly, the location of the point 
at which the boundary condition is satisfied is important 

According to Goldstein and Young (reference 6), 

11 in compressible flow the pressure increase at an;/ point 
of the body is l/(3 times the pressure increase in 
incompressible flow at the same point." That is, 

* ^ 

~ ( C M ) = — — ( e 0 ) 

U v ’ ' p u u ’ u 
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Comparison of this relation with equation (A6) shows that 
the Goldstein-Young method is also valid for two- dimensional 
oroblems but gives an incorrect result for three-dimensional 
problems. 
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